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Abstract 



We report on a study of the charged-skyrmion or spin-texture excitations 
which occur in quantum Hall ferromagnets near odd Landau level filling fac- 
tors. Particle-hole symmetry is used to relate the spin-quantum numbers 
of charged particle and hole excitations and neutral particle-hole pair excita- 
tions. Hartree-Fock theory is used to provide quantitative estimates of the en- 
ergies of these excitations and their dependence on Zeeman coupling strength, 
Landau level quantum numbers, and the thicknesses of the two-dimensional 
electron layers. For the case of near three we suggest the possibility of first 
order phase transitions with increasing Zeeman coupling strength from a many 
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skyrmion state to one with many maximally spin-polarized quasiparticles. 
74.60Ec;74.75.+t 
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I. INTRODUCTION 



In the limit of strong magnetic fields where all electrons are confined to a single orbital 
Landau level, interacting two-dimensional electron systems (2DES's) exhibit a rich variety 
of unusual properties. The quantum Hall effect, which occurs when the chemical potential of 
the 2DES has a discontinuity at a density n* which depends on magnetic field strengthB, is 
a prominent example. In the quantum Hall effect the incompressible groundstate at density 
n* can be a strong ferromagnet, i.e., its total spin quantum number S can equal N/2 so 
that electronic spins are completely aligned by infinitesimal Zeeman coupling. Surprisingly, 
in light of the strong magnetic fields, physics associated with spontaneous magnetization in 
these systems is experimentally accessible because the Zeeman coupling is typically quite 
weak compared to other characteristic energy scales and can even be tuned to zero, for exam- 
ple by the application of hydrostatic pressure to the host semiconductor. Recently there has 
been considerable interest in charged excitations of the incompressible groundstate at z/ = 1 
which has total spin quantum number S = N/2. As first noticed in numerical exact diago- 
nalization calculations,^ and demonstrated in recent experimentJiS the spin-polarization in 
these systems is strongly reduced away from u = 1 where the groundstate must incorporate 
the charged excitations of the = 1 state. This behavior can be understood,ili quanti- 
tatively in the limit of very weak Zeeman coupling, by identifying the elementary charged 
excitations with the topological solitons (skyrmions) of the 0(3) non-linear sigma (NLa) 
model in two spatial dimensionsB. With an appropriate kinetic termi, the latter describes 
the long wavelength T = dynamics of any Heisenberg ferromagnet in two dimensions. Two 
features distinguish the quantum Hall case. First, the skyrmions carry electrical charge as 
a consequence of their topological charge,Bi'0 and hence have a stable finite size for small 
but non-zero Zeeman coupling. Second, they are present in the groundstate near (but not 
precisely at) 1^=1, and as a consequence have an obvious influence on observable properties. 

In this article we discuss the elementary charged excitations of quantum Hall ferromag- 
nets using a Hartree-Fock approximation approach.i0 The Hartree-Fock approximation can 
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describe charged excitations which are sufficiently locahzed in space to invahdate the gradient 
expansion that underhes the NLa model description. A weakness of both the Hartree-Fock 
and NLa model calculations, which we will discuss further below, is the failure to respect 
the quantizationEUli of total electronic spin. The article is organized as follows. In Section 
II we make some general remarks on the implications of particle-hole symmetry for relation- 
ships between the energies and spin quantum numbers of positively and negatively charged 
excitations of quantum Hall ferromagnets and for the neutral decoupled particle-hole pair 
excitations which are important in activated transport experiments. El The considerations in 
this section do not depend on the Hartree-Fock approximation. In Section HI we explain 
some formal aspects of the Hartree-Fock approximation calculations we perform in order 
to explore the elementary charged excitations and how their energetic ordering depends on 
Zeeman coupling strength. Section IV presents and discusses our numerical results. We 
comment on the importance of the orbital Landau level index0 and on the thickness of two- 
dimensional electron layers. We conclude with a brief summary in Section V. An appendix 
discusses the spectral transfer in the presence of the skyrmion, analogous to that underlying 
midgap states in solitons in one-dimensional systems. 

We close this section with a remark on the terminology used in this paper. In NLa 
models the energy of a ferromagnet is expressed in terms of a function which specifies 
the direction of the local ordered moment as a function of the two-dimensional spatial 
coordinate. The energy functional in the "pure" NLa model invoked in discussions of broken 
symmetry states contains only a gradient term and the topological solitons of this model can 
be determined exactly.0 As the gradient term is scale invariant in two dimensions, the energy 
of these solutions is independent of their size. The NLa model appropriate for quantum 
Hall ferromagnets with a small Zeeman coupling has two additional terms which compete: 
a Zeeman coupling which favors small skyrmions and a Coulomb interaction which favors 
large skyrmions; together these determine the size and energy of the skyrmion and also its 
precise proffie which differs from those of the pure NLa model solutions^. In a microscopic 
quantum treatment^'i a spin-multiplet of elementary charged excitations with total spin 
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S = N/2 — K exists for each non- negative integer K. At large values of K, the correlations 
of the quantum states are well described by the classical soliton solutions; presumably, a 
treatment of the fluctuations about them, along the lines of soliton quantization in other 
problemsB, would reproduce the exact states. However, frequently the lowest energy charged 
excitations occur at relatively small integer values of K where both the symmetry restoring 
quantum fluctuations and the neglected higher gradient terms are large and the field theoretic 
description is no longer accurate. While, strictly speaking, one might wish to reserve the 
term "skyrmion" for excitations that are well described by the NLa model solitons, in this 
paper we take the liberty of referring to all elementary charged excitations of quantum Hall 
ferromagnets as skyrmions. 



II. PARTICLE-HOLE SYMMETRY 

When the spin-degree of freedom is included, the particle-hole symmetry of the 
Hamiltonian0 for interacting electrons in the lowest Landau level relates electronic states 
at Landau level filling factors v and 2 — v. Here v = N/N^ where N is the number of 
electrons and = AB/^q = A/(27r£^) is the orbital degeneracy of the Landau level. {B is 
the magnetic field strength, $o = hc/e is the electronic magnetic flux quantum, and i is the 
magnetic length.) The particle-hole symmetry of this system occurs because, apart from the 
constant quantized kinetic energy which is conventionally chosen as the zero of energy, the 
Hamiltonian contains only the term describing interactions of electrons within a degenerate 
Landau level. In the occupation number representation, many-particle states can be spec- 
ified either by the set of single-particle states within the Landau level which are occupied 
or by specifying the set which are empty, i.e. by specifying the states occupied by holes 
in the Landau levels. It is convenient to combine the particle-hole transformation with a 
spin-reversal so that the creation operators for spin-up particles are mapped to annihilation 
operators for spin-down holes and creation operators for spin-down particles are mapped to 
annihilation operators for spin-up holes. Under this mapping 
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N^^N^-N^^ N[. (1) 

It follows that N = N^ + Ni-^N' = 2N4,-N and S, = {N^ - Ni)/2 ^ S'^ = S, where S, 
is the z component of the total spin. The Hamiltonian H changes by a constant 

H = H' + 2{N~ N^)eo = H' + 2{N^ - N')e^ (2) 

where H' is normal ordered in terms of hole creation and annihilation operators and is 
identical to H except for the replacement of electron operators by hole operators, and eo is the 
energy per electron in the groundstate at = 1 which is readily calculated for any specified 
model electron-electron interaction. For 2D electron systems with a uniform neutralizing 
positive background charge 

^o = 4/(0^^P(-|^l'^V2M^) (3) 

where v{q) is the Fourier transform of the effective electron-electron interaction. 

In the following sections we will evaluate the energy change in the system when a single 
electron is removed from the system at a fixed magnetic field. As mentioned above and 
shown in Reference 0, the set of elementary charged excitations at = A^^ — 1 is composed 
of a single spin-multiplet with S = N/2 — K for each non- negative integer K. With finite 
Zeeman coupling, the lowest energy state in each multiplet has Sz = S = N/2 ~ K = 
N^/2 — {K +1/2) and the energy of this state relative to the u = 1 groundstate of the 
quantum Hall ferromagnet may be written as 

€K = UK + g*f^BBiK + 1/2). (4) 

Here g* is the g-factor of the host semiconductor, and Uk may be interpreted as the in- 
teraction contribution to the internal energy of the quasihole. It follows from the particle- 
hole transformation (Eq. (^) that the corresponding elementary charged excitations at 
N = N^ + 1 have S, = S = N' /2 - K = N^/2 - {K + 1/2) and energy relative to the 
N = N^ groundstate given by 
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et: = UK + 2e^ + g*^iBB{K + 1/2). (5) 

For example, the state with = is the single-Slater determinant elementary excitation 
which appears in the standard Hartree-Fock approximation. Note that and differ by a 
constant so that their minima always occur at the same value of K. For two-dimensional elec- 
trons in the lowest Landau level with purely Coulombic interactions, Uk=o = (e^/e£)(7r/2)^/^, 
and eo = -{e^ / e(!.){-n /^f^ so that e^^^ = Uk=o+9*I^bB/2 and e^^Q = g* 1138/2. It follow J 
from the NLa model that Uk^oo = (3/4)[/i^=o so that large K states will have lower energy 
if the Zeeman coupling is sufficiently weak. At intermediate values of K numerical calcu- 
lations of one sort or another are necessary to estimate Uk- In the following section we 
use a generalized Hartree-Fock approximation to obtain realistic estimates of Uk including 
effects of the finite thickness of the quasi two-dimensional electron layers and to discuss the 
influence of the Landau level index for which the quantum Hall ferromagnet occurs. 

The value of K at which and are minimized has been determined experimentally 
both by measuring the groundstate spin-polarizationiS near p = 1 and by measurin g0 the 
dependence on Zeeman coupling strength of energy required for the creation of free-particle 
hole pairs aX u = 1. For u close to one, interactions between elementary charged excitations 
can be neglected^l so that the partial spin polarization is given for z/ < 1 by 



and for z/ > 1 by 
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C^^^l-2{K + l){u-l)/u. (7) 



Note that d^/dv = 2K just below v = 1 and —2{K + 1) just above u = 1. The polarization 
^ may be measured in experimentsi'S that are sensitive to the spin magnetization of the 
system. Finally, the activation gap A measured in transport studiesEl is the energy to make 
an unbound particle-hole pair: 

A = e^ + e], = 2UK + 2eo + g*fXBB{2K + l). (8) 



The derivative of this energy with respect to the strength {g*fiBB) of the Zeeman couphng 
term is 27^^+ 1. As we discuss at greater length below, experiment^H'El performed at typical 
fields in GaAs 2DES's are all consistent with K = 3, in good agreement with calculations 
presented in subsequent sections. 



III. QUASIPARTICLE ENERGIES AND EXCITATION GAPS 

In this section some formal aspects of the Hartree-Fock theory for the the charged exci- 
tations at u < 1 are briefly discussed. Further technical details are provided in a previous 
paper on this topic! and in the appendix. (The excitations at z/ > 1 may be obtained 
from these by particle- hole symmetry as explained above.) Because of the symmetry of 
the skyrmion charged excitations, it is convenient to work in the symmetric gauge where 
single-particle wavefunctions in the lowest Landau level (LLL) have the simple form 

^ .-exp(-|.|V4£^) 
^""^ ' (2-+%£2^!)V2 • 

Here m = 0, 1, . . . , is the angular momentum, and z = x + iy expresses the 2D coordinate 
as a complex number. Note that states with larger m are localized further and further from 
the origin. We consider single-Slater determinant states of the form 

oo 

1^) = n (^-«ln + ^n^bl+l)\0 > (lO) 
m=0 

where |0) is the particle vacuum and |umP + |'i^mP = 1 so that the wavefunction is normalized. 
Here a}^ creates a down-spin electron and b}^ creates an up-spin electron in the mth angular 
momentum state. 

The form of these wavefunctions is essentially dictated by the symmetry of the classical 
skyrmion solutions which are invariant under the action of ± Sz for the skyrmion (an- 
tiskyrmion). This plus the requirement of LLL occupancy uniquely picks Eq. (|T0[) as the 
appropriate Hartree-Fock wavefunctions. It is also easy to demonstrate using Eq. (^) that 
the expectation value of the total spin-operator in this state describes a spin-texture with 
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unit topological charge provided that Um varies slowly with m from Um=o = 1 to Wm^oo = 0. 
{vm=o = and fm->oo = !)• 

This wavefunction is a generalization of the single-Slater determinant wavefunction 
proposed^ by Moon et. al. for which the expectation value of the vector spin operator 
gives a spin-texture identical with that of the pure NLa model skyrmion. The additional 
variational freedom in the wavefunction of Eq. ([10|) allows deviations from the pure NLa 
model so that microscopically localized charged excitations may be optimized in a way which 
depends on the details of the 2DES being considered. Far from the origin, this state is lo- 
cally identical to the ferromagnetic groundstates and all spins are aligned with the Zeeman 
magnetic field which is assumed to point in the 'up' direction. Near the origin the pro- 
jection of the total spin along the field direction becomes negative. It is easy to establish 
that the total decrease in electron charge near the origin compared to the ferromagnetic 
groundstate (11^=0 1 0)) corresponds to one electron. The total number of reversed spins 
in this wavefunction is 

oo 

K=Y.M'. (11) 

m=0 

The K = Hartree-Fock hole excitation is obtained by choosing Um = for all m. Note 
that the quantization of the number of reversed spins is not captured by these generalized 
Hartree-Fock variational wavefunctions. The quantization is obviously of some importance 
in the usual experimental circumstance since K is small. In this case better variational 
wavefunctions could be obtained, in principle, by projectingllllli the Hartree-Fock wavefunc- 
tions onto states with definite numbers of reversed spins. Here we account for quantization 
when necessary simply by restricting ourselves to solutions where K (the mean number of 
reversed spins) is an integer^. 

We find optimized skyrmion wavefunctions by minimizing (^|if|\I'), where the Hamil- 
tonian includes electron-electron interactions and Zeeman coupling. This procedure gives 
rise to a set of self-consistent equations listed in Ref. ^ which need to be solved numerically. 
Some typical results for the spin and charge densities relative to the groundstate are illus- 
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trated in Fig. |l], for skyrmions in the lowest Landau level in a sample of vanishingly small 
thickness, for several values o g, where g = g*^iBB is the Zeeman coupling strength!!. 
The total energy of this state, relative to the energy of the ferromagnetic groundstate, can 
be separated into interaction energy and Zeeman energy contributions as in Eq. (^): 

e-{K{g)) = U{K{g)) + ~g[K + \]- (12) 

For each value of g the self- consistent Hartree-Fock equations determine the value of K which 
minimizes e and, given K, the shape of the skyrmion state (specified by the Um values) which 
minimizes the interaction energy. For pedagogical purposes and for comparison with the field 
theoretical approach it is useful to consider the interaction energy as a function of K rather 
than g. U{K), which is just the Legendre transform of e~{K), can be obtained from the 
K{g) and e^{K{g)) produced by the self-consistent Hartree-Fock calculations using Eq. (|12|). 
Formally minimizing Eq. (|12D with respect to K we find that the optimal value of K for a 
given g is determined by 

dU{K) 



dK 



(13) 



Thus the global relationship between the equilibrium value of K and the Zeeman coupling 
strength is conveniently summarized in a plot of ^ versus K. We will use such plots in the 
next section to discuss the importance of the finite thickness of two-dimensional electron 
layers and of the Landau level index of the ferromagnetic state in determining K{g). The 
naive expectation is that the lowest energy skyrmion state should monotonically shrink 
as g increases; we see from Eq. (|13D that this is possible only if d'^U{K)/d^K is positive 
definite which is the usual requirement of convexity needed to obtain a continuous Legendre 
transform. The maximum value of —dU{K)/dK gives the maximum value of g at which 
Hartree-Fock skyrmion states with K occur. 



IV. HARTREE-FOCK THEORY NUMERICAL RESULTS 
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A. Lowest Landau Level: n=0 



In Fig. 1^ and Fig. ^we show results for e{K{g)), A{K{g)) and K{g) obtained by solving 
the Hartree-Fock equations at a series of g values for the case of a strictly 2D electron system 
and for the case of a quasi 2D electron system with the finite layer thickness modeled to 
approximate the typical experimental situation. Finite thickness can be introduced into 
the calculation by assuming all the electrons are in the some confined state x{^) of either 
the inversion layer or quantum well in which they reside. The effective two-dimensional 
electron-electron interaction may then be written in the form 

which is softer than the Coulomb interaction at short distances, but at large distances 
approaches the 1/r form. The precise wavefunction chosen for our calculation is the Fang- 
HowardU form x{z) = [UW^f^^^ ze''/^^ , where W may be understood as a measure of 
the thickness of the two-dimensional layer. Finite thickness corrections tend to reduce the 
energy scales of this system; for example, eo = —0.6267 for the strictly 2D model and 
eo = —0.3963 for the quasi 2D finite-thickness model for W = 0A5i. However, we have 
found that the optimal values for the total number of flipped spins in the skyrmion states 
are not dramatically affected by finite thickness. 

In Fig. 1^, we see that both e and A rise rapidly toward their K = values as g/{e^/etj 
increases. The rapid increase in the energy with increasing g is associated with the rapid 
shrinking of the optimal size of the skyrmion charged excitation, as seen in Fig. ^. For typical 
experimental systemsi'llll g ~ 0.015 — 0.020e^/e£ and the measured number of flipped spins 
i^' ^ 3, in good agreement with the present numerical results. The small Zeeman coupling 
strength required to reduce skyrmions to microscopic size is at first sight surprising and 
reflects the relatively weak K dependence of U{K). As explained in the previous section, 
U{K = 0) — U{K = oo) = (7r/32)^/^(e^/£) for the Coulomb interaction model. For large 
K the principal correction to the NLcx model K oo value for U{K) is the Coulomb 
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self- interaction of the excess charge. Assuming that the spin and charge structure of the 
skyrmion have the same size this energy is ~ [e^ j t^K'^l"^ . From Eq.( p!3| ) the optimal value 
of K implied by this approximation is oc g~'^^^ while the change in e is proportional to g^^"^ . 
A more careful analysis^ by Sondhi et al. replaces g in these estimates by g\ln(g)\. The g^^^ 
behavior at small g means that even weak Zeeman coupling compared to the characteristic 
interaction energy is sufficient to eliminate most of the energy difference between the lowest 
energy charged excitation and the K = Hartree-Fock quasiparticle energy. 

As explained in the previous section we can extract results for U{K), the internal energy 
of a skyrmion hole as a function of the number of reversed spins, from Fig. ^ and Fig. |^. 
These results are shown in Fig. ^ Note that (PU{K)/dK'^ is positive definite. The property 
that K{g) is not markedly altered by finite- width corrections to the effective interaction is 
associated with the weak K dependence of the difference between the two curves in this 
figure. The results shown here hint at a practical technical difficulty with these Hartree- 
Fock calculations. In practice we are forced to truncate the set of single-particle angular 
momenta we include in our wavefunctions for both the ferromagnetic groundstate and for 
the charged excitation at a finite value of m = Mmax- Since the single-particle orbital with 
angular momentum m is localized near a ring with radius £(2(m + 1))^^^) this corresponds 
to working with a finite size electron disk of radius R ^ £(2Mmax)^/^. The energy difference 
between the ferromagnetic groundstate and the state with the charged excitation will be 
given accurately by our calculation if the tail of the disturbance associated with the charged 
excitation does not extend to the edge of our system. As the Zeeman coupling weakens 
and the skyrmions get larger, finite-size effects become increasingly important. It is for this 
reason that the U{K) results shown in Fig. ^ do not accurately approach the analytically 
known K oo asymptote which for the case of strictly 2D interactions has the value 
(97r/32)^/^(e^/e^) ~ 0.94(e^/e^). (It is possible however, to get an accurate value for this 
limit by regularizing the ^ = skyrmion by placing it on a sphere. This is sketched in the 
appendix.) The smaller the value of g the larger the value of Mmax required to get accurate 
results. The convergence properties of our calculations with respect to Mmax are illustrated 
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in Fig. ^ where K is plotted as a function of {g\ log(5')|) From the asymptotic 




we know that for M, 



max 



oo this curve will follow a roughly straight line for large K. 



For each finite M^ax the numerical results bend away from this line, underestimating the 
number of reversed spins in the lowest energy charged excitation. From Fig. ^ we estimate 
that calculations performed with Mmax = 120 are accurate for g > 0.01 (e^/£), calculations 
with Mjnax = 480 are accurate for g > 0.005(e^/£) and that calculations with M^ax = 1920 
are accurate for g > 0.002(e^/£). The number of angular momentum states required to 
obtain good convergence for the total number of fiipped spins K is much larger than K, 
despite the fact that the Um and Vm approach their groundstate values relatively rapidly 
for m > K. As expected from the variational nature of the Hartree-Fock calculations, the 
energies presented in Fig. |^ and in Ref. approach their asymptotic values much more 
rapidly with increasing Mmax than estimates of the optimal value of K. Finally, we note 
that, except where the contrary is explicitly stated, results reported here and in Ref. ^ were 
obtained with Mmax = 120. In particular this means that the values of K presented in Fig. 
2 below g = 0.01 (e^/£) are slightly underestimated; however, we have found no significant 
errors in the energies of Fig. 1 for ^ > 0.002(e^/£). 

In Fig. ^ it is interesting to note that U{K) approaches K = with a finite slope. This 
is seen most clearly in Fig. ^ which, using Eq.( |T3|), can be regarded as a plot of K vs. 
~dU{K)/dK. The maximum slope of U{K) occurs at = and since \dU/dK\ increases 
monotonically with decreasing K this slope specifies the largest value of g for which a if 7^ 
solution of the Hartree-Fock equations occurs. The property that the Hartree-Fock skyrmion 
decreases continuously to zero size with increasing g contrasts with properties which would 
follow from other forms for U{K). For example if U{K) approached its i^' = value from 
below quadratically, the maximum value of —dU{K)/d{K) would occur at a finite value of 
K and the skyrmion would suddenly collapse to zero size once g exceeded this value. If 
U{K) approached its if = value from below as K"^ with s < 1, solutions with finite K 
would exist at arbitrarily large g. Finally if U{K) was an increasing function of K at small 
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reaching a maximum aX K = K*, solutions with K < K* would not exist at any g. We 
will see below that this is precisely the situation which often occurs when the quantum Hall 
ferromagnet groundstate has electrons in n 7^ Landau levels polarized. 

As remarked in the introduction, it is important to note that our Hartree-Fock approach 
does not respect the separate quantization of spin and orbital angular momentum in our 
system, and instead permits only the quantization of the difference of these two quantities. 
In the Hartree-Fock approximation, as in classical field theories, skyrmions appear as bro- 
ken symmetry states of the Hamiltonian. There is a similarity here to the standard BCS 
treatment of superconductivity, in which particle number is not a good quantum number in 
the mean-field groundstate0lli. However, while in the BCS problem there really is a broken 
symmetry in the thermodynamic limit, in our problem the skyrmion states only break the 
symmetry over a finite volume. Hence quantum fluctuations around the mean-field state 
will, if treated exactly, restore the individual quantization of spin and angular momentum. 
It is therefore important to address the extent to which these fluctuations influence the re- 
sults in Fig. 1^. In an exact calculation Fig. |^ would take the form of a series of straight line 
segments, since the quantum number K can only take on integer values^!. We expect the 
Hartree-Fock approach to yield a smooth interpolation of the exact results which is most 
accurate when its K value is integral. Comparison with exact finite size calculationJ^lifl 
confirm this expectation. (We note the fortunate feature that the Hartree-Fock treatment 
is exact both at vanishing g, where the infinite skyrmion is a classical object up to trivial 
global rotations and at large g where it describes the fully polarized quasihole.) 



B. First Landau Level: n = 1 



Quantum Hall ferromagnets occur at u = 1 and also at larger odd integral filling factors. 
In the Hartree-Fock approximation, the groundstate for u = 2n + 1 has all orbitals of both 
spins with Landau level index less than n occupied and only majority spin orbitals occupied 
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in the n-th Landau level. If Landau level mixing is ignored, the fully occupied Landau 
levels play no role and the theory of quantum Hall ferromagnets is altered only by the 
change in the index of the Landau level onto which the electronic Hilbert space is projected. 
The Hartree-Fock single-Slater-determinant state is an exact eigenstate of the Hamiltonian 
just as in the u = 1 case and is expected to be the groundstate, at least if n is not too 
large. The change of Landau level index may be accounted for without approximation 
simply by including an additional form factor correction to the effective electron-electron 
interaction. The Fourier transform of the interaction for electrons projected into the nth 
Landau level Vn{q) = v{q)[Ln{\)]^, where L„ is the nth Laguerre polynomial, and 
v{q) is the unprojected interaction. As for u = 1, the K = and K oo limits of the 
charged excitation energy can0 be calculated analytically when finite-thickess corrections 
are neglected. Wu and SondhiS have recently pointed out that for n > 1, Uk ^oo > Uk=o 
so that Hartree-Fock quasiparticles have lower energy than large spin-texture quasiparticles 
even when the Zeeman energy is not included. Recent transport experimentJii seem to 
imply that for u = 3, the lowest energy charged excitations have i^' = 0, in contrast with 
the n = case and in agreement with theory. In this section we examine the influence of 
the finite width of the electron layer on the energetics of the quasiparticles. 

The competition between Hartree-Fock and charged spin-texture quasiparticles can be 
understood in terms of the analytically known expressions for the Hartree-Fock exchange 
energy per electron for a full Landau level (eo in Eq.( ^) and for the spin- stiffnesffl 
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IQn J (27r)2 

Since Uk=o = — 2eo and Uk^oo = — eo + 47rps the difference is 

AU ^ U{K = 0) - U{K = oo) = -eo - 47rp, = 1 1 ^[2 - q']v^{q)e-^'/' . (15) 

For the Coulomb interaction modelS AU = {-l/16){TT/2y/^{e^ /i) in the = 1 (z/ = 3) 
Landau level and has a larger negative value for the n = 2 {u = 5) case. The [Li(g^)]^ form 
factor for electrons in the n = 1 Landau level strengthens the effective interaction at large q 

15 



which has more importance in ps than in eo because of the factor in its integrand. When 
finite thickness form factors are also included in the effective interaction, the importance 
of large q in these integrals will diminish. For sufficiently wide quantum wells it is clear 
that the sign of Af/ will be positive and charge spin textures will also occur in the n = 1 
Landau level. The objective of the calculations reported on below is to quantify the stability 
properties of skyrmions in n = 1 quantum Hall ferromagnets with finite width electron 
layers. We find that skyrmions in the first Landau level can be stabilized by finite thickness 
corrections, although for realistic widths the maximum g for which they are the lowest energy 
quasiparticles is roughly an order of magnitude smaller than for the n = case. Thus the 
observation of skyrmions in the first Landau level would require specialized samples or 
experimental techniques to access this very low Zeeman energy limit. 

Figure || illustrates the energy for both the K = Q quasiparticle (dotted line) and the 
optimal i^' > skyrmion (solid line) as a function of g for a typical experimental sample 
width, W = 0.45^. The skyrmion excitation is lower in energy for g < 2.4 x 10~^e^/e£. 
The results in Figure |^ suggest that for filling factors close to, but slightly away from, 
z/ = 3, a strong first-order phase transition should take place as g is increased from zero, 
in which the spin polarization changes very abruptly. That this transition should involve 
a very large number of flipped spins can be seen by noting that the optimal values of K 
of the skyrmions for small values g is extremely large, as illustrated in Fig. |^. Unlike the 
n = case, skyrmion solutions of the Hartree-Fock equations become unstable at a finite 
value of i^" (~ 9 for W = 0A5i) at which g = —dU{K)/dK reaches its maximum value 
(~ 0.003e^/e£ for W = 0.45£). For g larger than this value there are no K ^ solutions of 
the Hartree-Fock equations. 

We estimate the number of flipped spins in the skyrmion state for W = 0.45£ at the 
critical value of g to be approximately K = 14, using our Hartree-Fock approach with a 
system size of Mmax = 480. It should be noted that at finite skyrmion concentrations, espe- 
cially for such large skyrmions, a calculationEl including interactions among the skyrmions 
would give a more reliable estimate for the jump in the spin polarization at the transition. 
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We expect that such interactions would most likely reduce somewhat the magnitude of this 
jump. (Related finite size studies have been carried out by Jain and Wulll.) 

In general, the stability of > skyrmions is most easily assessed by considering the 
internal energy Uk- As in the previous section, Uk may be computed numerically by finding 

using the Hartree-Fock method, and then subtracting off the Zeeman contribution to the 
energy, as in Eq. ^ The results of such a calculation are illustrated in Fig. |^ for several 
values of the layer thickness W . These results confirm that skyrmions are the lowest energy 
charged excitations for large enough layer thickness, but become higher in energy than the 
K = quasiparticle as the system approaches the two-dimensional limit. For n=l, Uk may 
have local minima both at = and dX K ^ oo. For this reason it is possible for the 
Hartree-Fock equations to have two separate solutions for the same values of as shown in 
Fig. ^. (In the NLa model calculation^, the energy of the skyrmion branch increases from 
that of the ^ = solution in precisely the same fashion as in the LLL case for there is no 
distinction between the direct Coulomb interaction in the two Landau levels for skyrmions of 
divergent size. From the independent calculation of the lower energy polarized quasiparticle, 
it follows that the skyrmion energy must be non-monotonic with size and, in the absence 
of intervening additional extrema, must have minima at both microscopic and infinite size.) 
It is interesting to note that, unlike the case of skyrmions in the lowest Landau level, the 
curvature of Uk changes with increasing K (cf. Fig. H). As mentioned in the discussion 
following Eq. |13|, we do not expect to find stable skyrmion solutions for the convex regions 
of the graphs in Fig. ^. Indeed, in order to obtain values for Uk in the small K limit, it is 
necessary to add a fictitious term to the Hamiltonian of the form = <y{Sz — 5'°)^, where 
Sz is the operator for the total spin angular momentum of the state. Since this term couples 
only to the total spin of the system, it favors states with {Sz) ~ 5*°, but does not affect 
the optimal shape of the spin texture within the subspace of states with this expectation 
value for the spin. Thus by varying both a and 5° we can obtain states with small values 
of K, and their intrinsic energy may be obtained by subtracting the expectation value (Ha) 
from the expectation value of the perturbed total Hamiltonian. We note that this procedure 
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produces internal energies that join smoothly onto those obtained for larger values of K 
where the skyrmion states are stable for large enough layer thickness. 

Our results for this section are summarized in Fig. ^, by a "phase diagram" in g and 
W which shows where the lowest energy charged excitations are K = {] quasiparticles and 
where they are K ^ Q skyrmion charged spin-textures. Even for the most favorable layer 
thickness, W ~ 1.2£, the maximum value of g for which skyrmions are stable is outside the 
accessible range in many typical experimental systems. Nevertheless, we emphasize that if a 
system with very small effective values of g could be fabricated, for values of v close to (but 
not exactly equal to) 3 this line will separate states of maximal polarization {K = 0) from 
states with highly degraded polarizations due to the presence of skyrmion spin textures. 

V. SUMMARY 

In this work we have studied the charged skyrmion excitations of quantum Hall ferromag- 
nets within the Hartree-Fock approximation. Exact relationships between the energies and 
spin quantum numbers of quasihole and quasiparticle excitations energies near u = 1 were 
derived using particle-hole symmetry. An internal energy function U (K) was introduced to 
describe the relative stability competition between skyrmions states with different numbers 
of reversed spins, K. Locally stable skyrmion states with K flipped spins can occur for any 
Zeeman coupling strength only if (PU{K)/dK'^ > 0. In the presence of Zeeman coupling, a 
stable skyrmion states must satisfy dU{K)/dK = —g*^BB. Results for the dependence of 
the K value of the lowest energy skyrmion excitations on Zeeman coupling strengths were 
presented both with and without finite thickness corrections. Finite thickness corrections 
were shown to stabilize skyrmions in the n = 1 Landau level, where for zero thickness K = 
quasiparticles are lowest in energy at any Zeeman coupling strength. We have proposed that 
a first order phase transition between states with relatively large skyrmions and states with 
K = quasiparticles will occur at = 3 at a critical value of the Zeeman coupling g and be 
accompanied by a jump in the spin-magnetization of the electron system. 
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APPENDIX A: MIDGAP STATES AND THE INFINITE SKYRMION 

The most intuitive account of the physics of the skyrmions relies upon the spin Berry 
phases generated by adiabatic motion in a textured (effective) magnetic field arising from 
the magnetic exchange among the electrons themselves. Geometrical considerations show 
that the Berry phases due to a unit topological charge precisely mimic those of an additional 
quantum of orbital flux. The incompressibility of the quantum Hall fluid then implies that 
an additional electron (or fraction thereof for fractional quantum Hall states) must then be 
present in the region of the texture which is therefore a (dressed) quasielectron. 

A more detailed description of this is to note that "up" and "down" spins, i.e. the 
spins that respectively point parallel and anti-parallel to the local direction of the field 
in the texture, see oppositely signed fluxes. Consequently while the "up" spin (effective) 
Landau level gains one state the "down" Landau level loses one, i.e. the presence of the 
skyrmion causes the transfer of one state beween the upper (empty) and lower (filled) bands 
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at z/ = 1 thus allowing the extra electron to be accomodated in a low energy state, albeit 
at the energetic cost of texturing the spins. This is strongly reminiscent of midgap states in 
one-dimensional systems, i.e. polyacetylene@, where a background soliton causes states to 
appear in a gap and again there is the cost of creating the soliton which has to be overcome. 

Since this description is very much of the one-electron variety, it is instructive to see how 
it works out in the the Hartree-Fock treatment. The situation is clearest for the "infinite" 
skyrmion where the texture is arbitrarily slowly varying. It is very useful to regularize this 
limit by putting the system on a sphere where it takes the form of a purely radial texture, 
n{6, (f)) = f . Again we may deduce the Hartree-Fock description by noting that the classical 
texture is invariant under the action of J = L + S. In the LLL approximation this uniquely 
picks the wavefunction and fixes the eigenoperators of the Hartree-Fock Hamiltonian. For a 
sphere with 2q flux quanta threading the surface, the latter are: 

h- 
h- 

where 

Uj = 

In this basis the Hamiltonian and wavefunction are, 

Hhf= e-c]-Cj_- e+c]+Cj+ 

lil<g-i/2 \j\<q+l/2 

\^)hf= n (A3) 

lil<9+l/2 

where |0) is the vacuum state. The energies e± depend upon the details of the interaction 
and not just on the symmetries of the state. (In the limit of infinite system size these become 
degenerate with the corresponding up and down spin eigenvalues at (exactly) u = 1 as the 
skyrmion becomes locally indistinguishable from the ferromagnetic state.) The form of Hhf 
clearly shows that transfer of one state across the gap. Evidently, all we have used in this is 
that the states of Hhf must be classified only by their eigenvalues under J and hence must 
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^j4-i/2 + IjI < 9 + 1/2 

Vja]_^/2 - ^j^i+i/2' \j\ < ? - 1/2 



(Al) 



lq + l/2 + j , g + l/2-j 

and Vj = \ 



2q + l 



2q + l 



(A2) 



occur in the multiplets j + 1/2 and j — 1/2. This will therefore also survive the perturbative 
inclusion of Landau level mixing. 

For finite skyrmions, the picture is (unfortunately!) less elegant. The excess charge is 
now localized in a finite region and hence affects the Hartree-Fock eigenvalues in the core of 



the skyrmion. This can be seen in from Fig. |TU|, where we have plotted the eigenvalues, now 
on the plane, for skyrmions at two different values of g/{e^ /et). While the general scenario 
of a transferred state between two sets of levels still holds, the individual sets themselves 
become degenerate only in the large skyrmion limit. 

Finally, we note that putting the infinite skyrmion on the sphere allows us to compute 
its energy very accurately and thus circumvent the problems at small g noted in the text. 
The only trick here is that we need to subtract the Hartree self-interaction of the skyrmion 
charge, e^j{2^Jq) (y^g is the radius of the sphere), which is substantial for small system 
sizes. This greatly improves convergence and the known infinite system result is recovered 
to within a few percent already at system sizes with q ^ 10. 
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FIGURES 

FIG. 1. Charge and spin densities for quasihole skyrmions in the lowest Landau level for a 

system of zero thickness, evahiated for different vahies of g. Maximum angular momentum state 
kept in these computations was Mmax = 120 (see text). 

FIG. 2. Minimum hole creation energy (solid line) and transport activation energy (dashed line) 
in a v = 1 quantum Hall ferromagnet as a function of g. All energies are in units of (e^/el) and 
results are shown for both strictly 2D and finite-thickness models. Maximum angular momentum 
state in these calculations was Mmax = 120. 

FIG. 3. Number of reversed spins per hole in a v = 1 quantum Hall ferromagnet as a func- 
tion of g/{e^/e£). Results are shown for both strictly 2D and finite-thickness models. Note that 
for self-consistent solutions of the Hartree-Fock equations g = —dU{K)/dK. Maximum angular 
momentum state in these calculations was Mmax = 120. 

FIG. 4. Internal energy of a skyrmion hole J7 as a function of the number of reversed spins K. 
Results are shown for both strictly 2D and finite-thickness models. Maximum angular momentum 
state in these calculations was Mmax = 120. 

FIG. 5. Number of reversed spins in a skyrmion hole at weak Zeeman coupling as a function 
of finite system size. 

FIG. 6. Energy of the K = Q quasiparticle (dotted line) and the HF i^T > skyrmion (solid 
line) as a function of ^ for = 3 (n = 1) and a 2DES with thickness parameter W = 0.45£. The 
system size used for this calculation was Mmax = 480. All energies are in units of in units of /ei. 

FIG. 7. Number of flipped spins K of the HF if > skyrmion as a function of g, for 
= 3 (77, = 1) and a 2DES with thickness W = 0AM. Note that for self-consistent solutions 
of the Hartree-Fock equations g = —dU{K)/dK. The system size used for this calculation was 
Mmax = 480. All energies are in units of |e^. 
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Vertical dashed line indicates the value of g above which the skyrmion is unstable against 
the untextured quasihole. 

FIG. 8. Internal energy Uk vs. K for skyrmions in the n = 1 Landau level systems of 
varying thicknesses W. Skyrmions are the lowest energy quasiparticles for sufficiently weak Zeeman 
coupling if Uk-^oo < Uk=o- We find two distinct solutions to the Hartree Fock equations for 
W = 0.1£ in the vicinity K = 13, leading to a cusp in Uk; for larger values of W only one 
solution is found and the cusp becomes a smooth minimum. Where (fU{K)/dK^ < solutions to 
the Hartree-Fock equations can be found only by adding unphysical terms to the Hamiltonian as 
discussed in the text. The system size used for these calculations was Mmax = 480. Energies are 
in units of e^/ei. 

FIG. 9. Phase diagram for the dependence of the globally stable Hartree-Fock quasiparticle 
at = 3 (n = 1) on system parameters g and layer width W. {g is in units of e^/ei and W is in 
units of i.) The phase boundary separates parameter values for which the K = quasiparticle is 
lowest in energy from parameter values for which K > skyrmions are lowest in energy. Systems 
in which these parameters can be adjusted so as to cross the phase boundary should exhibit a jump 
in the magnetization. These calculations were performed with Mmax = 480. 

FIG. 10. The Hartree-Fock eigenvalues for skyrmions in the LLL at g/{e'^/e£) = 0.1 (stars) 
and 0.010 (circles). The former is the K = quasiparticle and the latter has K ^ 6. In both cases 
there is an extra state in the lower branch. As there are 301 orbitals in the system Mmax = 300, 
this means that the upper branches have 300 and the lower (occupied) branches have 302 states. 
Note that the gap is smaller near the origin which is where the extra charge resides and near the 
boundary due to loss of exchange. The shift in the eigenvalues near the origin is due to the Hartree 
repulsion of the extra charge. Evidently, the eigenvalues are much more uniform for the skyrmion 
but the overall reduction of the gap in its case is a finite size effect. 
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